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ABSTRACT 

An extended Hubbard model on a two-leg ladder is numerically studied by means of 
the quantum Monte Carlo techniques. The model we study has the nearest-neighbor 
interactions which are repulsive along chains and attractive for rungs. The plots of 
the doping parameter versus the chemical potential show two cliff-like regions and a 
large plateau region. Results on the charge susceptibility suggest it diverges in these 
cliff-like regions. These observations might imply a signal to the phase separation, 
which should be related to the effective attractive interactions along chains. 
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§1 Introduction 

The experimental realizations of the ladder systems|l|], which is considered as a 
useful bridge to clarify the relation between one- dimensional and two-dimensional 
systems, have given theorists exciting topics in the study of the strongly correlated 
electrons. Works on the quantum spin systems]^, to which the half-filled Hubbard 
system reduces in the limit of the strong repulsive Coulomb interactions, show a 
striking contrast between the ladders with even legs, for which the finite spin gaps 
exist, and the odd-legged ladders which, as well as the one-dimensional (chain) case, 
show no spin gaps. The t-J model on ladders is another object to be studiedQ in 
order to understand the doped system with the strong on-site Coulomb repulsion. 
Although we have not come to a very definite conclusion yet, interesting features 
like the ci-wave correlations are suggested. Furthermore some numerical studies for 
the Hubbard ladder model have been carried out|Q, but the issue of the long-range 
order of the pair-pair correlation functions seems controversial. 

Recently various extended Hubbard models are under study to find the effects 
of the longer range interactions between the electrons!^. In this paper we report 
numerical results obtained by the quantum Monte Carlo techniques on an extended 
Hubbard model with the two-leg ladder, whose Hamiltonian includes the Coulomb- 
like nearest-neighbor interactions that are repulsive along chains and attractive for 
rungs, studied in the conditions that the system is doped by holes or by electrons 
and the on-site Coulomb repulsion lies in the intermediate region. It should be noted 
that, for the values of the model's parameters we use in the simulation, the double 
occupancy is abundant so that holes can emerge even when electrons are doped to 
the system. Our motivation to study this model is to investigate behaviors caused 
by double-occupancy-pairs or vacancy-pairs on rungs, where (and hereafter) we refer 
to the state of a vacant site, a site upon which no electrons are sitting, as vacancy. 

We measure some thermodynamical quantities concerning to the number of par- 
ticles of the system such as the doping parameter and the charge susceptibility as 
a function of the chemical potential. We observe that the plots of the doping pa- 
rameter versus the chemical potential show two unexpected cliff-like regions at finite 
dopings in addition to a few plateau regions, which have been expected from the re- 
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suits for the standard Hubbard model on the chain or the square lattice systems[|]. 
The charge susceptibility measured in each of the cliff-like regions shows a large 
peak which seemingly diverges at the temperature consistent with zero to imply the 
possibility of the phase separation that has not been confirmed in the preceding 
studies of the Hubbard-like models. 

In §2 we make a brief description on the model and method. The Monte Carlo 
results are in §3 together with the definitions of the quantities we measured. Finally 
54 is devoted to discussions. 



§2 the Model 

The Hamiltonian we study is 

Nr-l Nr 

cr I i=l a i=\ 

Nr-\ Nr Nr 

a, a' I i=l a, a' i=l I i=l 

where Ci^i^a denotes the annihilation operator for an electron with spin a (t or |) 
which is located on the i-th rung along the leg / (a or b) of a ladder, A^^ being the 
total number of rungs and nj_; g. = c| ;^^Cj_;^o-- We employ an open boundary condition 
to avoid the hoppings between the sites on the first rung and the sites on the A/^^-th 
rung. 

The partition function Z is given by 

Z = tr{e-'3(^-'^)}, (2) 
with the inverse temperature (3, the chemical potential fi and 

Nr 

■A^^EEE^M,., (3) 

o- I i=l 

which is the total number of electrons sitting on the ladder. 

In the Suzuki- Trotter formula[0 with the finite Trotter number N^, 

^ ^ E E • ■ ■ E(«i I e-^^/"^* I «2) ■ ■ ■ {ajv. I e-^^/^' I (4) 

ai Q2 aivj 
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where we denote 7i — fiAf by J'. We employ a complete set 0] 

I a) =1 5*1, 5*2,- --SNr) (5) 

with sixteen states in Table I to denote the state Sj on the i-th rung, which are the 
eigenstates for 

o" cr,a' I 

^ <t,<t' I 

so that the checkerboard decomposition is available for J = J^-, 

Nr/2 Nr/2 Nr 

Jl = -tcY^HY. [4fc-l,«,aC2fc,Z,<x + h.c] +VcY,Y.Y. n2k-l,l,an2k,l,a' + X^^^' 
o- I k=l a,cT' I k=l i=l 

Nr/2-1 Nr/2-1 

J2 = J2 [4k,l,aC2k+l,l,a + h.c] + K X! H J2 ^2fc,«,a^2fc+l,«,a' 

o- / k=l a,a-' I k=l 

Nr 

+ Ej'^, (7) 

i=l 

where an even Nr is assumed. Parameters ui and U2 in Table I are given by 
1 



Ml = 



U 1 

1 + , U2 = 



§3 Numerical results 

Now let us show our quantum Monte Carlo results on the model described in 
the previous section. Here we concentrate our attention on particular values of 
parameters in the Hamiltonian (P, = 1, = 2, Vc = 2, V,. = —4 and U = 4. 
For this choice it is possible to obtain statistically meaningful results while the sign 
problem is quite serious for other choices we tried. 

Because of the emergence of negative weights, we should measure a physical 
quantity A in each thermally equilibrated configuration together with its sign and 
obtain its expectation value by 



where Z^{Z_) is the number of configurations with positive (negative) weight, and 
y4+(A_) denotes the contribution from positively (negatively) signed configurations. 
The r ratio defined by 

r = — , (9) 

will inform us how serious the cancellation is. 

Changing the number of rungs Nr {8 < Nr < 32), the inverse temperature f3 
(3.3 < P < 10) and the chemical potential fi, we measure the doping parameter 

(10) 

where = (A/')av is the number of electrons and A^^ = 2A^r. the number of sites on 
the ladder, and the charge susceptibility per site. 

Typically we generate 80000 configurations in one measurement by the local update 
and the global update in the Trotter direction, details of which are in 0, and last 
60000 configurations are used for the measurement. Each datum is an average of 
six such measurements with different random number sequences and the statistical 
error of the datum is calculated from the standard deviation among these six runs. 
As we observed very little discrepancy between the Nf = 32 results and the Nt = 64 
ones we fixed A^t=64 throughout the measurements. The r-ratio in the simulations 
varies from 0.17 to 0.78, depending on values of Nr, fi and /?. 

Figure 1 presents the doping parameter 6 on Ng = 16 (namely A^^ = 8) and 
Ng = 64 ladders at /5 = 10 as a function of the chemical potential fi. We observe 
rapid increases of 6 around ~ — 1, where the system is under the hole doping, 
and fi ^ 5 which belongs to the region of the electron doping, together with three 
plateaus in between. Comparing the results for A^^ = 16 and Ns = 64 ladders, we 
conclude that these plateaus would be split due to the finite size effect and would 
merge to form the 6 = plateau in the thermodynamic limit, implying a large 
charge gap Ac ~ 6. It seems that 6 changes smoothly in the neighborhood of the 
half-filling. 
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Let us then present the data on Xc the charge susceptibihty per site, in more 
detail. Figure 2(a) plots Xc per site versus 5 at /3 = 10 on up to an Ng — 64 ladder, 
choosing several values of ji with which we observe two outstanding enhancements 
of Xc corresponding to the results in Fig. 1. It is notable that for each value of /i 
in the negative 5 region the values of Xc and 5 strongly depend on the size of the 
ladder, while the data for positive 5 are much less size-dependent except for the 
value of Xc around 5 ~ 0.2. In both the positive and the negative 5 regions the 
thermodynamic limit of the peak value of Xc for each value of /3, which we denote 
by xT-: is successfully obtained by the hnear extrapolation of the data 1/xc at the 
peak for each ladder size versus l/Ng with Ng = 16, 24, 32, 48 and 64. In Fig. 2(b) 
we plot these extrapolated values as a function of the temperature T = l/j3. (Note 
that the ordinate measures the inverse of xT-) Linearly fitted lines to the data 
drawn in the figure suggests that l/xT reduces to zero (namely xT S^^^ to infinity) 
near T = in both regions of 6. Figure 2(c) shows the values of 6 to maximize Xc 
evaluated from the data in Fig. 2(a), which we denote by 5o, as a function of the 
1/Ns at temperatures T = 0.1, 0.2 and 0.3. Because of large uncertainties compared 
to those on the Xc at peaks it is difficult to estimate reliable critical value Sc at which 
the charge susceptibility diverges, but it is intriguing to note the data suggest finite 
critical values (5c ~ 0.2 and 5c ^ —0.2. 

§4 Discussions 

In the previous section we showed our quantum Monte Carlo results on the ex- 
tended Hubbard ladder model described in §2. What we measure in the simulation 
are the quantities related to the particle density, the doping parameter 6 and the 
charge susceptibility per site Xc which could signal the charge gap and the phase 
separation. We find two regions where Xc is greatly enhanced, one of which corre- 
sponds to the electron doping {S < 0) and the other to the hole doping {5 > 0). 
Linear extrapolations on these enhancements suggest that Xc diverges for the dop- 
ings away from the half-filling at the temperatures which are consistent with zero. 
Between these singularities we observe three plateaus which would merge to one in 
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the thermodynamic hmit to indicate the existence of the large charge gap. 

We observe no sign for the singularity of Xc near 5 = 0, which is in marked 
contrast to the observations in the one- and two-dimensional standard Hubbard cases 
where Xc behaves like 1/5 when 5 ~ 0[^ |TU|. Since it seems difficult to understand 
our results in the same manner as the latter cases, this feature, together with the 
enhancement of Xc at finite values of 5, implies that a new picture is necessary to 
explain the behaviors of our model. We would like to stress that this is not quite 
surprising because studies of the spin 1/2 ladder systems, which corresponds to the 
U ^ oo limit of the half-filled Hubbard models on a ladder, present a picture that 
differs much from those for the spin 1/2 models on a chain or on a plane [|l], 

What happens to the model in these regions with the increasing Xc around the 
non-zero 6c then? Although it is an open question so far and much more work would 
be necessary to answer it, one clue might be found in the behavior of the double- 
occupancy-pairs on rungs in the positive 6 region as well as the vacancy-pairs on 
rungs in the negative S region. 

Let us comment on the hole doping case first. Here we examined the number of 
doubly occupied sites A^^ and the number of double-occupancy-pairs on rungs Ndd 
defined by 

Nr 

iVd = (A/'d)av, AAd = 51 E^'M.T^'Ma (12) 

l=a,b i=l 

Nr 

Ndd={J^dd)a.v, J^dd = ^ni^a^ni^a,['^i^b,^ni^b,[^ (13) 

i=l 

respectively. In the measurement of these quantities for the hole-doped A^^ = 64 
ladders at /3 = 10, we observe that the ratio of doubly occupied sites to all sites, 
Nd/Ns, linearly increases from to ~0.4 as 6 decreases from ~ 0.5 to ~ 0. Compar- 
ing this ratio with the number of double-occupancy-pairs on rungs divided by the 
number of rungs, Ndd/Nr, we see that most of these doubly occupied sites prefer 
pairing on rungs. In order to clarify contributions of the double-occupancy-pairs 
on rungs to the charge susceptibility Xc, we then divide Xc into three parts Pi, P2 
and P3 and measured them separately. Here we let Pi represent the portion which 
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is proportional to the susceptibility concerning to the double-occupancy-pairs on 
rungs, Xdd, 

Pi = - {^Afdd)l) = X 8xdd, (14) 

Xdd^^a^dd')..-{Afdd)l), (15) 
P2 the portion concerning to other states, 

P2 = j^iiKth').. - {Kth)l), Kth = M- aUm. (16) 

and P3 the interference term 

P3 = 2j^{{MddNoth)^. - {Mdd)^.{-^oth)e..)- 

Figure 3 shows results for Pi and P2 on an A^^ = 64 ladder with /3 = 10 together 
with clS cl function of the chemical potential /x. We clearly see that the double- 
occupancy-pairs on rungs have dominant effects in the rapid increase of Xc- 

In the case of the electron doping we turn our attention to the number of the 
vacant sites, A*"^, calculated by 

K = N,-Nd- {N, - 2Nd) = NJ + iVrf, (17) 

the number of the vacancy-pairs on rungs, N^y, and the susceptibility concerning to 
these pairs, Xw, which are defined using /ij ^ o- = 1 —rii^i^^ instead of nj^o- in eqs.(|l^) 



and (15). We observe approximately particle- hole symmetric behaviors of A^^ and 



Ny^ for —0.5 < 6 < —0.1 in comparison with the hole doping case. We also observe 
that the susceptibility Xm shows a remarkable enhancement which suggests that 
the vacancy-pairs play an important role when Xc diverges in the region 6 ~ —0.3, 
although the relation between Xw and Xc is more complicated than that between 
Xdd and Xc- 

Several remarks on the future work are in order. Correlation functions corre- 
sponding to the susceptibilities Xc, Xdd and Xw would be worth studying to under- 
stand whether the long-range order exists in this model. It will be important to 
measure the magnetic susceptibility and the specific heat in order to shed the light 



8 



on the physics of the model. As was mentioned in §2, we had to hmit ourselves to a 
particular choice of parameters tc, tr, Vc, K and U in the simulation because of the 
sign problem. It would be a future task to study this model for other values of the 

parameters. It would be also interesting to investigate such an extended model in a 
bilayer system whose planes are the two-dimensional counterparts of the chains in 
the ladder system, where the realistic choice of parameters might be realized. 
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Figure Captions 
Figure 1 

Monte Carlo results on the doping parameter 5 as a function of the chemical potential 
fi measured on A^^ = 16 and Ns = 64 ladders at the inverse temperature f3 = 10. 
All the statistical errors are within symbols. 
Figure 2 

(a) The charge susceptibility per site, Xc, versus the doping parameter 6 measured 
for the chemical potential ranging —1.35 < < —0.95 and 5.0 < < 5.3 on the 
Ns =16, 24, 32, 48 and 64 ladders at /? = 10. All the statistical errors for 6 are within 
symbols. The statistical errors for Xc are within symbols if not shown explicitly. 

(b) Temperature dependence of the inverse peak value of the charge susceptibility 
per site extrapolated to A^^ oo, which we denote by l/xT- Squares (circles) are 
symbols for the data with negative (positive) 6. Errors evaluated in the process of 
the extrapolation are within symbols. The result from the linear fit for the data 
with negative (positive) 6 is also shown in the figure with a solid (dotted) line. 

(c) Values of 6q, the values of 6 to maximize Xc at T = 0.1, 0.2 and 0.3 as a function 
of 1/A^s- We evaluate 6q by the peak among the measured data in Fig. 2(a) and 
estimate its errors by the difference of the two values of 6 around the peak. 
Figure 3 

Contributions of the double-occupancy-pairs on rungs and other states to the charge 
suseptibility Xc measured in the positive 6 region on an A^^ = 64 ladder at the inverse 
temperature jS = 10, as a function of the chemical potential fi. Open circles, open 
squares and filled circles represent Pi, P2 and Xc defined by eqs. (0),(|l6]) and (p4]), 
respectively. Statistical errors are within symbols. 
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Table Captions 
Table I 

States on two sites of the i-th rung used to construct a complete set in the Suzuki- 
Trotter formula. | 00) represents a state with no electron on either site of the rung. 



No. 




1 


1 00^ 


2 




3 




4 


1 fc^ , + 1 00^ 


5 


1 fc^ , - 1 00^ 


6 


r'^ tr^ t 1 00^ 


7 


;^(4a,T4M+4aa4T)|00) 


8 




9 




10 




11 




12 


Tf^'^U.T'^laa'^lfe,! 4,a,t4,6,t4,fe,J,) 1 00) 


13 


^i4,a,'[4,a,A,b,f ~ Ag^^aAaO I 00) 


14 


'^i4,a,t4.,a,l4.,b,l + 4,aa4,6,T4,6,i) 1 00) 


15 


'^i4,a,'l4,a,l4,b,l ~ 4,a,l4M4,b,i) 00) 


16 


4,a,A,a,A,b,A,b,l 1 00) 



Table I 
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